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COUPLING AND HARNACK INEQUALITIES 
FOR SIERPINSKI CARPETS 

Martin T. Barlow and Richard F. Bass 

Abstract. Uniform Harnack inequalities for harmonic functions on the pre- and 
graphical Sierpinski carpets are proved using a probabilistic coupling argument. Var- 
ious results follow from this, including the construction of Brownian motion on Sier- 
pinski carpets embedded in R'', d> 3, estimates on the fundamental solution of the 
heat equation, and Sobolev and Poincare inequalities. 



The Sierpinski carpets (SCs) we wiU study are generahzations of the Cantor set. 
Let Fo = [0, l]"* be the unit cube in R'', d > 2, centered at Zo = (1/2, . . . , 1/2). 
Let k, a be integers with 1 < a < fc and a + k even. Divide Fq into k"^ equal 
subcubes, remove a central block of a'' subcubes, and let Fi be what remains: 
thus Fi ~ Ff) — ((fc — a) /2k, {k + a)/2k)'^. Now repeat this operation on each of 
the fc'' — a"^ remaining subcubes to obtain F2. Iterating, we obtain a decreasing 
sequence of closed sets F„; then F = fX^^aFn is a Sierpinski carpet and has 
HausdorfT dimension df — df{F) — \og(k'^ — a'^)/log(fc). (When d = 2, fc = 3, and 
a = 1, we get the usual Sierpinski carpet.) Let Fn = k"Fn C [0,ooY, and define 
the pre- Sierpinski carpet hy F — U^i (see [10]). The graphical Sierpinski 
carpet is the graph G — (V, E) with vertex set V = (zq + Z"*) n F and edge set 
E={{x,y}€V:\x~y\^l}. 

Thus int(F) is a domain in R'' with a large-scale structure which mimics the 
small-scale structure of F. We are interested in the behavior of solutions of the 
Laplace and heat equations on F, F, and G. One reason for this is applications to 
"transport phenomena" in disordered media (see [6]); another is the new type of 
behavior of the heat kernel on these spaces. Let W be Brownian motion on F with 
normal reflection on dF, and let q{t, x, y) be the transition density of W, so that q 
solves the heat equation on F with Neumann boundary conditions on dF. 



1991 Mathematics Subject Classification. Primary 60B99; Secondary 60J35. 

Key words and phrases. Harnack inequality, Sierpinski carpets, self-similar, fractals, Brown- 
ian motion, heat equation, transition densities, Poincare inequality, Sobolev inequality, spectral 
dimension, electrical resistance. 

Received by the editors November 9, 1992 

Research partially supported by NSF grant DMS 91-00244 

©1993 American Mathematical Society 
0273-0979/93 $1.00 -I- $.25 per page 



1 



2 



M. T. BARLOW AND R. F. BASS 



Theorem 1. There exist ci, . . . , cg € (0, oo) and ds = ds{F) e (1, df) such that if 
x,y G F, t € (1, oo), |a: — y| < t, then 

cit-'^^/^ exp ( - C2 (J f—) j 

<q{t,x,y)<C3t-'^'^^expl-C4[^ j, 

where d^ = 2df/ds; while if x,y G F, t & (1, oo), |a; — j/] > t, then 

(2) exp(-C5 j <q{t,x,y) <exp{-ce 

The index dg is called the spectral dimension of F and turns out to be much more 
significant than the HausdorfF dimension dj as far as analytic properties of these 
spaces are concerned. Since dg < d, (1) confirms the physical intuition that the 
presence of increasingly large reflecting barriers causes heat to dissipate to infinity 
more slowly. It seems unlikely that there is any simple relationship between dg, k, a, 
and d. 

While there is a well-developed approach to the heat equation using analytic 
tools such as Sobolev or log-Sobolev inequalities (see [7]), these methods do not 
appear to give the best-possible results on spaces such as F — compare the upper 
bound on q{t,x,y) given in Theorem f with the results of [10]. 

The proof of Theorem 1 rests on the following Harnack inequality. Let D cM.'^ 
be open: we will say that h is harm,onic on D C] F ii (i) Aft = in int(D n F) 
and (ii) h has normal derivative a.e. on D n dF. Equivalently, h is harmonic with 
respect to W{t A To), where To = inf{i : W{t) ^ D}. Let D„ = (-1, fc")"^. 

Theorem 2. There exists ci e (0, oo) {depending only on d, k, a), such that if h is 
positive harmonic in Dn flF and x,y G Dn-i r\ F, then h{x)/h{y) < ci. 

Remarks. 1. Note that c\ is independent of n; otherwise the result is trivial. 

2. The case d = 2 was proved in [1]; the proof there relies on the fact that a 
closed curve in the plane separates the plane into two pieces. Just as in the case of 
elliptic operators, the results for two dimensions are considerably easier to prove. 
The result of [1] was extended in [8] to SCs with ds{F) < 2. 

3. Using the symmetry of F, Theorem 2 extends to other domains in F. 

4. Theorems 1 and 2 actually hold for a much wider class of SCs, those satisfying 
a higher-dimensional generalization of (2.1) of [4]. 

5. A similar result holds for the graphical Sierpinski carpet G. 

6. Most existing proofs of Harnack inequalities for selfadjoint operators depend 
on Sobolev inequalities, which in turn depend on the underlying geometry of the 
space. Here the appropriate Sobolev inequality involves the spectral dimension 
ds{F); however, no geometric definition of dg is known. Thus we were led to aban- 
don analytic approaches in favor of the probabilistic coupling argument described 
at the end of this paper. 

We now describe some other consequences of Theorem 2. Let J^ = U^o k'^F, the 
SC extended to [0, oo)"^, and write ^ for Hausdorff a;"^^ -measure on F. 
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Theorem 3. There exists a strong Markov process Xi with state space F such 
that X has a strong Feller transition semigroup Ft which is ^-symmetric, Xt has 
continuous paths, Xt is self-similar with respect to dilations of size fc", and the 
process X is locally invariant with respect to the local isometrics of F. 

Let p(t, X, y) be the transition density of Xt with respect to /i. Then p(t, x, y) is 
the fundamental solution to the heat equation on F: du/dt = Aj,u, where is 
the infinitesimal generator of Xt- Then we have 

Theorem 4. There exist ci, 02,03,04 e (0, 00) and dg = ds{F) G {^,df) such that 
for all x,y £ F, t & (0, 00), 

,<-./^exp(-c,(fcft)'"'"-') 

< ,(t, X, ») < c,t--'^ exp ( - c. ""'-""), 

where = 2df/ds. Moreover, p{t, x, y) is C°° in t, andp{t, x, y) and all its partial 
derivatives with respect to t are jointly Holder continuous in x and y. 

Many properties of the process X, such as its transience or recurrence, the 
existence of local times, the existence of self-intersections, and the asymptotic fre- 
quency of eigenvalues follow easily from Theorem 4. For example, note that X is 
point recurrent if and only if ds{F) < 2. 

The next set of consequences include Sobolev inequalities, Poincare inequalities, 
and electrical resistance inequalities for F, F, and G — nine theorems in total. Since 
the electrical resistance inequalities are probably the least well-known type, we give 
the one for G as a representative sample. If B is any subset of G, let \B\ denote 
the cardinality of B. Then R{B), the resistance from B to infinity, is defined by 

R{B)-' = inf J (•/■(^) - /(y))' : / = 1 on f{x) ^ as |a;| ^ 00 I . 

[{x,y}eE(G) J 

The inverse of R{B) is the conductance from B to infinity and equals the capacity 
of B. 

Theorem 5. Suppose dg = ds{F) > 2, and let ( = ds/{ds — 2). Then there exists 
ci such that if AcG, \A\ < cii?(A)~^. 

Theorem 5 follows fairly straightforwardly from Theorem 4 by applying ideas of 
[11] and [12]. As Theorems 1, 3, and 4 follow from Theorem 2 by generalizations 
and modifications of methods of [1-4, 8, 9], we discuss only Theorem 2. 

Let Wt be the Brownian motion on F described above, and let 

t{x, r) = mi{t : \Wt - x\ > r}, T{x, r) = mf{t : \ Wt - x\ < r}. 

The following lemma is proved in a similar fashion to Lemma 3.2 of [1]. 

Lemma 6. There exist C2 > Ci > 1, 5 > independent of r such that if x,y G F 
and \y — x\ < cir, then 

(3) P^(T(x',r) < r(x,C2r)) > (5. 

It is known (see Theorem 3.9 of [5], for example) that the Harnack inequality 
Theorem 2 follows from (3) and an oscillation inequality of the following form. 
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Lemma 7. There exists p < I such that if n > 1 and h is positive harmonic on 
Dn n F, then 

(4) \h{x) - h{y)\ < p sup Jh{z)\, x, t/ e £>„_! n F. 

zeDnCiF 

To show (4), it suffices to construct two F- valued Brownian motions and 
W^, starting from x and y respectively, which couple (i.e., meet) with probability 
at least 1 — p before either exits Dn- 

Fix n. Let Sm be the collection of cubes of side length fc™ with vertices in 

k"^!!^. Say that x, y E F arc m-associatcd if there is an isomctry of the cube in Sm 
containing x onto the cube in Sm containing y that maps x onto y. Note that if 
two points are m-associated, then they will also be ^-associated for all £ <m. 

Suppose first that x and y arc m-associated. We start a Brownian motion W^{t) 
onFatx. Let C/q = 0, and t/j+i = inf{f : |M/^(t)-H^^(f/i)| > fc™}. The key step is 
to exploit the local symmetry of F to construct, using suitable reflections, another 
Brownian motion VF^(i) on F, starting at y, such that (a) W'^{t) and W^^(f) are 
m-associatcd for alH > 0, and (b) there exist j and ci > such that 

(5) P(l^^(f/j(a;)) and Wy{Uj{uj)) are (m + l)-associated ) > ci. 

A renewal argument and then an induction show that if x and y are 0-associated, 
then W-^ and couple with probability C2 > before either process leaves 
Lemma 6 then follows easily. 
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